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Abstract. It is known that for any non-zero M G ( — |, oo), if we roll the conic {(x, y) : Ax^ ~ 1/ = 
1} on a line in a plane, and then we rotate about this line the trace of a focus, then we obtain a surface 
of revolution D(M) with mean curvature 1. If M < 0, D(Af) is embedded and it is called unduloid, 
if M > 0, D(M) is not embedded and it is called nodoid. These surfaces are called Delaunays and 
they are foliated by circles. The trace of the focus in the construction above is called the profile 
curve of the Delaunay and it is transversal to the circles. Another well known family of constant 
mean curvature surfaces are the Twizzlers, they are foliated by helices and we can naturally define a 
profile curve which is transversal to the helices. To make the presentation in this abstract easier, we 
will be only considering Twizzlers with mean curvature 1. In this paper we will prove that if we roll 
the profile curve of a Twizzler on a line in a plane and, simultaneously, we move the points in the 
line at the same speed, so that the rolling motion of the profile curve looks like if it were placed on 
a treadmill, then, the trace made by the origin of the profile curve is one of the closed heart-shaped 
curves ]HI§(M, w) = {{x, y) : x^ + + = M} for some M > - j and w > 0. We have that, 

up to a rigid motion in R"', every Twizzler is determined by its corresponding heart-shaped curve 
]HIS(M, w); we will denote this Twizzler by T(M, w). We will prove the continuity of the map p from 
fl — {M + vi £ C : M > — J, M ^ 0, 0<u<l}to the set of immersions that sends, every point 
of the form — ^ + to a cylinder of radius | ; a real point A-I to D(M) and, for any i\/ > — | and 

V > 0, the point M + vi to T(A/, \J^^)- This map p is one-to-one in the interior of Q,. Surfaces in 
the image under p of the purely imaginary points in Q., that is, Twizzlers associated with M = 0, 
contain the axis of symmetry and they are special in the sense that each one of these Twizzlers is 
isometric to a nodoid and to an undoloid, more precisely: the undoloid p{~-fj^^^), the twizzler 
p{v\) and the nodoid pijj^^W)^) isometric. Moreover, we prove that for any u G (0, 1), all the 
surfaces in the image of the curve a„(M) = p{M + \ ^J^^'^^IzImXIwI+Tm^^^ ) defined in the 
interval (— „ , ,-, ) are isometric. Notice that this curve starts at an undoloid, then it 

passes through an Twizzler that contains the axis of symmetry and then it ends with at a nodoid. 
We will prove that every Twizzler that contains the axis of symmetry is contained in the interior of 
a cylinder of radius 1. Moreover, it is either properly immersed with a discrete group of rotations 
acting on its isometry group, or it is dense. As an additional result, we prove that if we do the 
same rolling on a treadmill procedure to the profile curve of a flat surface with the same helicoidal 
symmetry as the Twizzlers, then, the trace of the origin must lie in a line. We finish the paper 
by showing the .pdf version generated by the software Mathematica of some of the programs that 
provide some pictures and animations showing property of Twizzlers and Delanunay surfaces. 



1. Introduction 

Let us start this section with some figures, most of them taken from animations created with 
Mathematica. The code that produced Figure 1.8 and the videos suggested in Figures 1.3, 1.6, 1.7, 
and 1.15 wih be included at the end of the paper. The reason why this code works is developed 
throughout the paper. 
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Courtesy of Ivan Sterling 



Figure 1.1. Dynamic interpretation of the profile curve of an unnodoid 




Courtesy of Bruce Solomon 
Figure 1.2. Dynamic interpretation of the profile curve of a nodoid 




Figure 1.5. The profile curve, Twizzler and the heart-shaped curve 




Figure 1.8. Fundamental piece of the profile curve 




Figure 1.10. Profile curves of properly immersed Twizzlers that do not contain the axis 




Figure 1.12. This graph shows how the maximum and minimum distance from the 
origin to the profile curve of T(M, v) changes with respect to M 
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Figure 1.14. Isometric Twizzlers 




Figure 1.15. Isometric nodoid, undoloid and special Twizzler 
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Figure 1.16. Points in the moduli space that represent isometric Twizzlers 





Figure 1.17. A surface with helicoidal symmetry is flat when its treadmillSled is a 
vertical semi line 

In 1841 [2], Delaunay showed that if one roll a conic section on a line in a plane and then rotates 
about that line the trace of a focus, one obtains a cmc surface of revolution. When the conic is 
a parabola we obtaine a Catenoid, when the conic is an ellipse, the surface is embedded and it is 
called an undoloid and when the conic is a hyperbola the surface is not embedded and it is called a 
nodoid. The undoloids and nodoids are called Delaunay surfaces. Figure 1.1 illustrates the relation 
between the ellipse and the trace of its focus. Notice that only one focus is used to get the curve that 
needs to get rotated in order to generate an undoloid. Figure 1.2 illustrates the relation between 
the hyperbola and the trace of the foci. Notice that both foci are used to get the curve that needs 
to get rotated in order to generate a nodoid. 

In 1970, Using the integrability of the Gauss Equation and Mainardi-Codazzi equation, Lawson 
showed that for any immersion /o : C/ — t- R"^ with constant mean curvature H, defined in a simple 
connected surface U, there exists a 27r-periodic 1-parametric family of immersions {fg : ?7 — t- : 
6 E [0,27r]} with constant mean curvature H and with the same induced metric. The family is 
called the associated family to /q. 

Remark 1.1. The map 6 — > fg is continuous with respect to the parameter 9. 

In 1982 [T] Do Carmo and Dajczer studied surfaces with constant mean curvature that are invariant 
under the group gt : R'^ — t- R^ of rigid motions 



5( (x, y, z) = {x cos(t) + y sin(t), —x sin(t) + y cos(t), z + ht) h G (— oo, oo) 

They called these surfaces helicoidal with a pitch h. When /i = 0, the group gt becomes a group 
of rotations and the helicoidal surfaces become surfaces of revolution. Therefore, when /i = the 
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helicoidal surfaces are either cylinders, spheres or Delaunay surfaces. The main result in their paper 
established that every Helicoidal surface is associated with a Delaunay surface. 

A Twizzler is an immersion of the form 



(1.1) (pi^s, t) = {x{s) cos{wt) + z{s) sin(u;t), t , —x{s) sm{'wt) + z{s) cos{wt)) 

with constant mean curvature. We will assume that the curve {x{s), z{s)) is parametrized by arc 
length and we will call it the profile curve of the Twizzler. Notice that Twizzlers correspond to those 
helicoidal surfaces with nonzero pitch. 

In this paper we will give an interpretation of the profile curve of Twizzlers similar to the interpre- 
tation for the profile curves of the Delaunays. In order to do this, given a regular curve a in the 
plane, we define the TreadmillSled of a as the trace of the origin when this curve (and the plane 
that contains it) rolls on a treadmill located at the origin which is aligned along the x-axis. Notice 
that the origin, as a point in the plane that contains the curve a, plays and important roll in this 
definition, for example, the TreadmillSled of a circle with center at the origin is just a point, while 
the TreadmillSled of a circle of radius R whose center is at a distance r of the origin, is a circle 
of radius r with center (0, i2). Figure 1.3 shows the TreadmillSled of an ellipse with center at the 
origin. The black dot represents the center of the ellipse while the blue curve represents the trace 
of this center. As the next theorem shows, it turns out that for any curve (x(s), z{s)) parametrized 
by arc-length there is an easy expression for its TreadmillSled. 

Theorem 1.2. If a{s) = {x{s)^z{s)) is a curve parametrized by arc-length and 



(,i{s) = x{s) x' (s) + z{s) z' (s) and (,2{s) = —x{s) z' (s) + z{s) x' {s) 
then the TreadmillSled of the curve a is — ('^i(s), ^2('S)) 

Notice that (^i(s), ^2('5)) are the coordinates of the vector (x(s), z{s)) with respect to the orthonor- 
mal basis {{x' (s) , z' (s)) , {— z' (s) , x' (s))} . 

There is a relation between the profile curve of Delaunay surfaces and conies. For Twizzlers, there is 
a relation between their profile curves and the level sets of the function h^{x, y) = x^^y^^—j^M^^. 

It is not difficult to check that the range of the function is the interval [— |, oo), that h^{^—^ = 
{(0, -i)} and that, every M > — j is a regular value of and h^^{M) is a closed simple curve. 
We will refer to these level sets as heart-shaped curves. Figure 1.4 shows some of these level sets 

The following Theorem establish the relation between the heart-shaped curves and the profile curve 
of the Twizzlers 

Theorem 1.3. The treadmillSled of the profile curve of a Twizzler with constant mean curvature 1 
is a heart-shaped curve —h~^{M) for some M > —j. The value M = is achieved by a cylinder 
of radius ^. 

The previous theorem is a consequence of the following proposition: 

Proposition 1.4. The immersions given by have mean curvature 1 if and only if the func- 

tions ^1 and ^2 defined in Lemma satisfy the following ordinary differential equations ^^(s) = 
/i(6(s),C2(s)) and ^'2(8) = f2{^i{s),C2{s)), where, 



-w'^X2 + 2(1 -I- w'^xl)2 
l + w'^{xl + xl) 



fl{xi,X2) 



X2 + 1 
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Moreover, the function is an integral curve for every solution {£,i{s) , £,2(3)) ■ 

Figure 1.5 shows the profile curve of a Twizzler and the heart-shaped curve associated with it. 
Figure 1.6, is a picture showing how, for this Twizzler, the TreadmillSled of the profile curve is 
indeed negative the heart-shaped curve. 

The following theorem gives us another application of the treadmillSled 

Theorem 1.5. A surfaces of the form is flat if and only if either the TreadmillSled of the 

profile curve is a point in the y-axis different from the origin, in this case the surface is a cylinder, 
or the treadmillSled of the profile curve is contained in a vertical semi line that starts at point in 
the X-axis different from the origin. Moreover, the functions x and z can be explicitly computed and 
they are: 



, , 1 2^/as'+b I . 2\J as + h , 2^/ as + b 1 . 2^/ as 

x{s) = -cos hvas + osm and z[s) = w as + cos sm ■ 



2 a a a 2 a 

Figure 1.17 shows how the TreadmillSled of the profile curve produces a vertical semi line. 

If we exclude the cylinder and the value M = — |, Theorem 11.31 establishes a 1:1 correspondence 
between pairs {M,w) with M > and vu > and Twizzlers with mean curvature 1. In order to 
deal with a bounded parameter, we will define v = + w'^), therefore, the parameter v moves 
from to 1 when w moves from 00 to 0. Figure 1.7 shows pictures from an animation that produces 
a piece of the Twizzler associated with values of M and v. We will refer to this Twizzler as 1^{M, v) 
when the dependence of M and v is needed. 

Explicit parametrizations for the heart-shaped curves make the understanding and drawing of Twiz- 
zlers easier. The following Lemma provides a possible parametrization for these heart-shaped curves. 

Lemma 1.6. For any M > — | and vu > 0, the curve a{t) = {pi{t), P2{t)) defined on the interval 
[0, 27r] and given by 



, , ^ . ^ , , , -1 + Jl + AM + Bcos^{u) sin('u) 

Pi(u] = Acosiu) and P2(u] = : 

^ ^ ^ ^ ^ ' 2^\^uflA^ cos2(n) 



where, 



-1 + Mw'^ + Vl + (1 + 2M)';i;2 + M2i(;4 

A = — — and 

V2w 



^ 2 + 2M'^w^ + w'^ + 2(Mu;2 - l)Jl + {I + 2M)w'^ + W^w^ 
B = 



is a closed simple regular curve that parametrizes the heart-shaped curve h^^{M). 
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We can define a fundamental piece of a profile curve as a connected piece of profile curve with the 
property that the TreadmillSled motion of this piece goes exactly once over the heart-shaped curve. 
It is not difficult to see that the whole profile curve is the union of fundamental pieces. Figure 1.8 
shows the fundamental piece of the profile of a properly immersed Twizzler, along with the whole 
profile curve made up of four pieces in this case, and the graph of the Twizzler. For an undoloid we 
could define a fundamental piece as the trace of the focus of the ellipse when this ellipse rolls once. 
It is clear that the integers Z acts on the group of symmetries of the Delaunay surfaces in the form 
of translations. Theorem 11.81 shows that the group Z also acts on the set of isometrics of Twizzlers. 

In order to produced Figure 1.8, it was necessary to find a formula for the length of the fundamental 
piece of the profile curve. The following lemma provides such a formula. 

Lemma 1.7. The length of the fundamental piece of the Twizzler 1^{M,v) is Jq^ ^J^du where, 

Ku) = fUpliu), P2(.U)) + f^{pi{u), P2{u)). 



The functions pi, p2, fi and /2 are defined in the Lemma \1.6\ and the Proposition \1.4\ Recall that 
w and V are related by the equation v = 1/(1 + w'^). 

The following theorem along with Theorem 11.121 and Theorem 11.131 give a precise picture of the 
Moduli space for the Twizzlers. 

Theorem 1.8. Every Twizzler %{M,v) is invariant under a group of rotations of the form {R{n9) : 
n G Z}, where the angle 9 depends on the fundamental piece of the profile curve. Moreover, if 
R[mO) = R{9) for some integer m, then the Twizzler is properly immersed, otherwise it is dense 
in the interior of a cylinder of radius 1 when M = or it is dense in the region bounded by two 
concentric cylinders of radii ri(M) = | Vi+4M- _i | r2(M) = in ifiQ case that M ^ 0. 

We also have another type of density: the set of points (M, v) associated with properly immersed 
Twizzler is non countable and dense. In addition, we have that T(M, v) is a properly immersed 
surface with profile curve consisting of b fundamental pieces, if and only if J^^ ipdu = 27r| with a 
and b positive relatively prime integers, where 



For a Twizzler that does not contain the axis of symmetry, the "properly immersed vs. dense" 
property established in the theorem above was proven by Hitt and Roussoss [1]. By numerically 
solving the equation ipdu = Itx ajb in the previous theorem, we can graph profile curves of 
Twizzlers with any desire property. In Figures 1.9 we solve the numerical equation fixing M = 
and taking several integer values for a and b. Since M = 0, these profile curves represent twizzlers 
that contain the axis of symmetry. In Figure 1.10 we took several values for M 7^ and a and 
b integers to produce properly immersed twizzlers that do not contain the axis of symmetry. In 
Figure 1.11 we took a and b so that 0/6 is not rational so that the twizzler is not properly immersed. 
In Figure 1.13 we took a = 5, 6 = 4 and 4 values of M in order to produce properly immersed 
examples, we also show the points (M, tt;) associated with these Twizzlers. Figure 1.12 shows the 
graph of the function ri(M) or r2(M) defined in Theorem 11.81 

By Dajczer and Do Carmo's Theorem, we know that every Twizzler is isometric to a Delaunay 
surfaces, in order to decide which one, we will prove the following properties for Delaunay surfaces. 
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Theorem 1.9. For every nonzero real number M E {—^,00), the Delaunay surface D(M) generated 

by the conic {{x,y) : — ^ = 1} has constant mean curvature 1. Moreover, the quotient between 
the maximum value of the Gauss curvature and the minimum value of the Gauss curvature o/B(M) 
is given by the function rs{M) = — (i^^^^S)^ 

We have that the function rs defines a bijection between the intervals (— j,0) to (0,1) and also 
it defines a bijection between the intervals (0, cxd) to (0,1). On the other hand, each undoloid is 
isometric (associated with angle ^ = ^) to a nodoid, see for example Appendix A in the paper [5] 
written by Kopouleas. As a consequence of Theorem 11.91 we have the following corollary: 

Corollary 1.10. Two helicoidal surfaces are isometric, if and only if the quotient of the maximum 
value and the minimum value of the Gauss curvature is the same, in particular, for any u € (0, 1) 
the undoloid D(— ,_, ) is isometric to the nodoid D( ,_, ^^,0 ) 

As pointed out before, each nodoid is isometric to a undoloid, therefore we can replace the word 
Delaunay by either the word undoloid or nodoid in the Do Carmo-Dajczer theorem, that is, we can 
say that each Twizzler is isometric (or associated) to either a nodoid or an undoloid. Another family 
of surfaces that holds the same property is the set of Twizzlers that contains the axis of symmetry, 
that is, the set of Twizzlers corresponding to M = in the moduli space. We will call these surfaces 
special twizzlers and we will denote them by §T(f), that is, ST(v) = T(0,v). 

The following Theorem gives us the quotient of the maximum value and the minimum value of the 
Gauss curvature for special Twizzlers. 

Theorem 1.11. For every nonzero real number v € (0, 1), the quotient between the maximum value 
of the Gauss curvature and the minimum value of the Gauss curvature of the special Twizzler surface 
§T(f), is —V. Moreover EZ{v) is isometric to the undoloid 9 (—-ij^^^j) ^'^d the nodoid J!)[ ^-^ J^^2 ) ■ 

Figure 1.15 shows two sets of isometric nodoid-undoloid-special Twizzler surfaces. We can generalize 
the theorem above as follows. 

Theorem 1.12. If v = jq^, the quotient between the maximum value of the Gauss curvature and 
the minimum value of the Gauss curvature of the Twizzler surface T(Af, v) is — . 
Moreover, fixing c G (0, 1), all the Twizzlers in the set 

fljg -'-'^'^ -(fljg + ^ + ^^) : M.(-^-^^,^-^)} 
^ ^ Vl + 4M + 1 - 2M + c(Vl + 4M - 1 + 2M) ' ^ {1 + {1 - ' ^ 

are isometric. 

Since the curve rv ( M) - (M v^l+4g-l-2M+c(v^l+4M+l+2Af ) x , ■ f, , ( ^5 s _ ( 

bmce tne curve ac{M} - [M, ^^-^^^_^^^^^^^^-^^^^^^^^) satishes tliat a^i (1+^)2 j - [ (1+^)2, 

ac{0) = (0, c) and ctdjj^^^) = (^j^^p-, 0) then, we have as a corollary of the previous two theo- 
rems and Remark 11.11 that 



Theorem 1.13. Let n = {M + xv £ C : M > -I, M ^ and < v < 1}. The function p from 
17 to the set of immersions in given by 

p{M + w) = T(M, v) for any v > and M ^ -\ 
p{M) = B(M) for any M ^0 and M ^ -\ 
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is continuous in the sense that for every point p in Q there exists a neighborhood U of p in Q 
and a continuous function f : U x — > such that for any M + G U, the map {s,t) — > 
f{M + iv, s, t) defines a parametrization of the surface p{M + if). Moreover, the function p is one 
to one in the interior ofQ.. 

The continuity at the points of form — ^+\v follows from Theorem ll.81 because each Twizzler T(M, v) 
is contained the region bounded by the two concentric cylinders of radii ri[M) = | vl+^^nij g^j^^j 

r2(M) = Figure (1.17) shows the trace of the curve au for several values of u. As a 

consequence of the previous theorem we have that the family of associated surfaces for Twizzlers 
repeats every ^. More precisely, it is vr-periodic and it is 2:1 when defined in the interval [0,7r]. 

The author would like to express his gratitude to Professors Robert Kusner, Ivan Sterling, Bruce 
Solomon, Wayne Rossman, Martin Kilian and loannis Roussos for solving several of his doubts 
about Delaunays and Twizzlers and for providing him references. 



2. Proof of the results 



In this section we will provide the proofs of the results established in the introduction. We will start 
with Theorem 11.21 

Proof Theorem II. 2t For every sq in the domain of the curve a, we can find an isometry from 
R^ to R^ of the form Tx = ii(7(so))x + c(so) that rotates the curve a so that a(so) moves to the 
origin and the vector a'(so) moves to the unit vector that points to the positive direction of the x 
axis. Here 



_ ( cos(t) -sin(t)\ 
^^'^ ~ Vsin(t) cos(t) ) 

The angle 7(so) and the constant c{sq) can be found by the solving the following equations (the 
solution of these equations is the key part of the proof of this theorem) 



/5'(so) = and /3(so) = 

where /3(s) = R{'^[sQ))a{s) + c{sq). If follows that the TreadmillSled of the curve a is the curve 
c(s). A direct verification shows that 




This completes the proof of the theorem. ■ 

Proof of Proposition II. 4t Since the curve [x{s),z{s)) is parametrized by arc-length, we can 
consider a function 6[s) such that 

x' {s) = cos(0(s)) and z {s) = sin(0(s)) 
Let us define the functions Ci(s) and ^2(5) as 
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(,1 = X cos{6) + z sm{6) and ^2 = —xsm{6) + zcos{6) 
A direct verification shows that 



(2.1) X = 6cos(0) -6sin(0), z = ^1 sin(0) + 6 cos(0) and 9' = x'z"-z'x" 
Moreover, it is not difficuh to check that, 

e'i = ^'6 + l, ^2 = -^'6 and f^+f^=x^ + z^ 
A direct verification shows that the first fundamental form of (p is given by 

E = {(j)„<l)s) = 1, F = {<l)sAt) = w{zx' -xz') = wC2, G = {<l)ucj)t) = l+w\x^+z^) = l+w\Cf+Ci) 
and therefore, 

EG-F^ = l + w\^f + el) - w^C^ = 1 + w^^f . 
A Gauss map of the immersion cj) is given hy u = (pg x (pt- A direct verification shows that, 

^{s, t) = {sm{wt — 6{s)),w(,i,cos{wt — 0{s))) 



Again, a direct verification shows that the second fundamental form of (p is given by 

e = ((/>,,, u) = — , f = ((/>of , u) = — , and q = (Sft, u) = — ^"^ - 

' ^/i+w^ei a/tt^ yrwef 

Therefore, if we assume that the mean curvature '^2{EG-F'^f ^ obtain the following ODE 



(2.2) 



-w^^2 + 2{l + w^^l)l 



Using the expression above for 9' in the equations = ^'^2 + 1> '^2 ~ ~^'?i) we obtain that .^i 
and ^2 satisfy the following ODE, 

(2.3) e'l = /l(6,e2) ^2 = f2{xi,X2) 

where, 



-w'^X2 + 2(1 + w'^xl)' 
1 + w'^{x'l + xl) 

w'^X2 - 2(1 + w'^xf) 
1 + u;2(xf + X2) 



jl{Xi,X2) = , , ..2^ 2 I 2:2 + 1 



w ^2 — ^y-L -r uj ^Yl' 



A direct verification shows that if we define 
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/i^ : R as hw{xi,X2) = == + xj + 

y 1 + w^xf 

then, hyu is a first integral of tlie ODE for and i-e. for any solution ^i(s) and (,2(8) of this 
system, we have that 



hw{(,i{s),^2{s)) = M where M is a constant 
This completes the proof of the Proposition. ■ 

Proof of Theorem II. 5t If we define the functions 9, .^1 and ^2 as in the previous theorem, we 
have that the equation: Gauss curvature equal zero, that is, the equation eg — = reduces to 
9' = Replacing this equation in the equations ^[ = 9' ^2 + 1, ^'2 ~ ~^'Ci) obtain that ^1 

and ^2 satisfy the following ODE, 



(2.4) = = ^ 

It follows that ^i(s) = ^ for some real number a. If o = then ^2 is also constant different from 
zero, and the surface (/> is a cylinder. In the case that a is not zero then ^2 = ±\/ai~-F6 and 
61(5) = =p2y^^±6^ 'j^j^jg completes the proof of the theorem. ■ 

As mentioned before. Theorem 11.31 follows directly from Proposition 11.41 Lemma 11.61 is a direct 
computation. 

Proof of Lemma 11.71 and Theorem II. 8t Let us start by pointing out that for any M > — i, 

the minimum and maximum distance from the origin to the set h~^{M) are ri(M) = | Vi+4Af^ | 
and r2(M) = respectively. We can see this by using the method of Lagrange multipliers 

to find the maximum value and minimum value of the function R{xi,X2) = + subject to 
the restriction hw = M. We can prove that if {xi,X2) is one of these critical points, then xi = 
otherwise the equations 



dR , dhw , dR , dh^, 

= — = — 

0x1 0x1 0x2 0x2 

has not solution. Once we know that xi must be zero, we obtain from the equation hw = M that 
X2 is either _ qj- \/^+iM~ ]_^ -^^^ continue with the rest of the proof. Let us assume that 

(x(s),y(s)) are such that the surface (jl.ip has constant mean curvature 1. Since the curve {pi,p2) 
defined in Lemma [l. 61 is regular, we have that 



au au 

is a periodic positive function. Likewise, since the functions /i(a;i,X2) and /2(xi,X2) only vanish 
simultaneously at (xi,X2) = (0, — |), we get. 



is a positive periodic function. Notice that ^i(s) = and ^2(5) = — ^ is the only constant solution 
of the system (12. 3p . For any other solution, since is a first integral of the system, there exists 
M > and a function a{s) such that 
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= Pi{cr{s)) and ^2(5) = P2(o-(s)) 
is a solution of the system (I2.3P . From the equations above we have 



(2.5) eiisf + e2isf = Xia{s))a'{sf 

On the other hand we have that, 

e'i(s) = fiiCi{s),Us)) = fi{pi{Hs)),p2{s)) and C'^{s) = f 2(^1(8) , Us)) = f2{pi{a{s)), p2{s)) 

It foUows that fj is either strictly increasing or strictly decreasing, WLOG we can assume that a is 
strictly increasing. Therefore we get that. 



a'{s) 



IpHs)) 



! A(a(.)) 

If n{u) is the inverse of the function (t{s), we have that 



(2.6) k'{u) 



'A(n) 



(t'(k(u)) y /i(u) 

If we change from the variable s to the variable u, that is, if we consider the functions, 

0(u) = 6{k{u)) il{u) = il{n{u)) i2{u) = Ul^iu)) X{u) = xiK{u)) Z{u) = ziK{u)) 

It follows from Equation (I2.6P and Equation (I2.2p . that 



9'(n)=^(n) where ^(,) = -^^(n) + 2(1 + u.V?(n))^ /A(u) 



i + w^{pI{u) + pI{u)) y p{ 



u 



Since ip{u) is a periodic function with period 2tt, it follows by the existence and uniqueness theorem 
of ODEs that if e{27r) = 60, then 

(2.7) 6{u + 2j7r) = JOq + 9{u) for any integer j 

Since \ {x{s), z{s))\ = |(^i(s),.^2('S))|, we have that the piece of profile curve 

^fP = C fundamental piece = {{x{u), z{u)) : U G [0, 27r]} 

also satisfies that ri(M) = min{| q\ ■ q G Cfp} and r2(M) = min{| q\ ■ q £ Cfp}. Using ()2.ip and 
(prp we get that 



(2.8) I f !! ! S:') = Hi ( fir!) where - " 



The equality above implies that the image of the profile curve can be viewed as the orbit of the 



group {Rg acting on Cjp, that is. 
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(2.9) C = {{x{t),z{t)) : i G R} = {R^^V : J e ^ and p G C fundamental piece} 

It follows from the equation above that li ^ is a. rational number, then C is a properly immersed 

curve and, if |^ is irrational, then C is dense in the annulus {(xi, X2) : ri{M) < y^x'f + x'2 < r2(M)} 
when M 7^ 0, or it is dense in the circle of radius 1 when M = 0. Therefore, Twizzlers with constant 
mean curvature 1 have the following property: they are properly immersed, or they are dense in the 
region contained between two concentric cylinders or they are dense in the interior of a cylinder of 
radius 1. We can prove that a surface corresponding to an irrational value |^ is dense by showing 
that the profile curve is dense, and we can prove that the profile curve is dense by showing that the 
intersection of this curve with a circle centered at the origin is either the empty set or dense in the 
circle. The problem of proving this last statement reduces to that of showing that for any irrational 
number l the set {l — [nt] : n £ 7^} is dense in the interval [0, 1] which is a known fact. To finish 
this proof we notice that since the function (x{s), z{s)) is parametrize by arc-length, we have that 

the length of the fundamental piece is k{2tt) = f^^ ^J~^^du. Also, since 0'{u) = Tp{u) we have that 
^0 = Jq'^ ilj{u)du. ■ 

Proof of Theorem II. 9t Let us assume that ID)(M) is parametrized as 



(j){s,t) = {x{s), z{x) sm{t), z{s) cos (t)) 

where the profile curve {x{s), z{s)) is parametrized by arc-length. A direct verification shows that 
if 9{s) is a continuous function such that x'{s) = cos{9{s)) and z'{s) = sm{9{s)), then the mean 
curvature of D(M) is ^{9' — — 7^^)- Since the mean curvature of D(M) is 1, we have the functions 
6{s) and z{s) satisfy the following ordinary differential equation 

e' = 2 + ^^ and z' = sm{9) 

This ODE has as a first integral the function h{z, 6) = z{cos{9) + z). Recall that the function z{s) is 
always positive. Since the minimum of the function /i is — |, it follows that there exists a non-zero 
constant k > — | such that h{z{s), 9{s)) = k. When k < 0, the level sets of h{z, 9) are bounded and 
therefore D(M) will represent an undoloid, when A; > the level sets are not bounded and D(M) 
will represent a nodoid. In any case, the z-values of the level sets of 9) are bounded. A direct 
computation shows that the maximum and minimum of the z-values of the level set h{z^ 9) = k are 
i+V^+4fc g^j^j I i--Vi+4fc I ^ prove that k must be equal to M by comparing these critical values 

of z{s) with the maximum and the minimum value of the profile curve viewed as the trace of the 
focus of a conic when it is rolled on a line. A direct computation shows that the Gauss curvature 
is _^_£^W g^j^j since 9' = 2 + £2^^ then the Gauss curvature reduces to _ '^°'^(^)(^^+'^°^(^)) ^ Using 
the Lagrange multiplier method we get that the maximum and the minimun of the Gauss curvature 
subject to the constrain h{z, 9) = k are 

4Vl + 4/c , Ay/l + Ak 

and 



(l + Vmk)^ (-1 + Vl + 4A:)2 

respectively. It follows the quotient between the maximum of the Gauss curvature, and the minimum 
of the Gauss curvature is — (^^ii^^S)^. Since k = M, the theorem follows. ■ 

Proof of Theorem I1.12t Using the same notation as in the proof of Proposition 11.41 we get that 
the Gauss curvature K satisfies 
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EG-F^ {i + w^il? + + 

Taking piiu) and p2{u) as in Lemma [1.6 1 we get the following expression for the Gauss curvature in 
terms of the parameter u 

-Avo^ Vl + 4M + B cos2 u sin u 



4 + + AA^w^ cos'^ u - 2w'^Vl + AM + B cos^ u sin u + (1 + AM)w'^ sir? u + w'^{SA^ + Ssin^ u) cos^ ■ 

A direct computation shows that the derivative of the function K = K{u) is of the form cos(n) po{u) 
where po{u) is a positive function, therefore the maximum of the Gauss curvature occurs when 

u = ^ and it is equal to „ ^I". X i and the minimum of the Gauss curvature occurs when 
u = and it is equal to — 2_|-(i^^2A/^ vTRm)^^ ' therefore we get that the quotient of the maximum 

value of the Gauss curvature and the minimum value of the Gauss curvature is — ^^f!"*"^^! ^i-^^^r -? • 

2+(l+2A/+Vl+4M) 

The last expression in terms of v transforms into 

1 + 2M - ^/l + AM + v- 2Mv + v y/l + AM 
1 + 2M + VI + 4M + v- 2Mv - v ^/TTAM 



A direct verification shows the last expression reduces to c when we replace v by ^iill^ \ ^^{"'"'^^^llLljl"'"^"'"^''^^^ . 

^ 1- J y'i+4M+l-2A/+c(Vl+4A/--l+2Af) ' 

therefore, all the Twizzlers T(M, ^/j-ip||j^j"2M+c(VTT^ ^ isometric. This concludes the 

proof. ■ 

Theorem 11.111 follows from Theorem 11.121 by taking M = 0. 
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3. Some programs using Mathematica 

In this section we will present the file with extension pdf generated by the software Mathematica 
for the programs that produce some of the figures in this paper. 

Mllllt= NTMtegi'ate [S^l; [25 Sin[t] '"2 • > Cub [t] - 2] , ft, 0, 2 7r}] 
OutllDl= 25.5:7 

-. a = 2/ »; 

b = 2/ Z5; s = NBSa3.VE [{at ' [t] == Cos [tli[t] ] , i ' [t] == Siii[tk[t] ] , 

th ■ [t] == (dC«s [th[t] J • 2 + b Siii[th[t] ] - 2) / (a K[t] Sin[th[t] ] -b z[t] Cos [th[t] ] ) , 
>c[ai = 0, z[0] - th[0] (X, E, tk} , (t , 0, 25. S2e»$SSC3102S3S }] r 

a = -{Ji[t] , z[t] )j 

It = { {Cob [tk[t] ] , Sin[th[t] ] } , { - Sin[tk[t] ] , Cos [ti[t ] ] }} ; 
f = Evaluate [ - R. a J . s] } 

pf = FarawtricFlot [f , {t, 0, ») , PlotlUiige -> Kll] ; 

tm = ParaneteicPlnt [{H-2, - .5) + .45 {Cos [t] , Sin[t] {2. - .5) + . 45 {Cos [»] , S±n[t] ) ) , 

({-2, - .5} + . 5 {C«H Hi. + 2 7t) / 1] , SlBL(t + 2lt)/2J}, {1, -.5} + .5 {Cm [t y 2] , Sin [t. / 2]} , 
{2ty7r, .5-. 5}, {it. / w, - .5-. 5)), {0, .45 - .5} + . 02 {Cos {t] , Sla{t.] ) ) , 
{t, -TT, 7T) , RxeB -F EiLLse , FlotStyle -t { lliifi , Directive [R«d, Thich] , Black}] r 
Hajilpiilate [Shav{ {FdEAmetrxcPlot {Evaluate [ 

E^atten {K / .s/.t-fi,l].i]i + Elatten {-It.a/.a/.t-*i, 1] I . s] , 
(t, 0, 25. 52«333tC3TB2595' ), FiotKange -+ Ull , 

PlotStyle -r Directive [Ked, Thickness [0.**!] ]] , pf , Graphics [Pout [Flaitten { f /. t -> 1] ] ] , 
tm, Graphics []lrfoir[{ { .S, - .42-. 5), - . 42-. 5J )]])], {1, », 25. 52C33tt( 3702595' ) ] 

Clear [ 
a, 

SI,;. 



Figure 3.1. This code produces the video suggested in Figure (1.3) 
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b2 = NDSoIv« [|lt ' [tj Cos [t*[t] ] , 

z ' [t] == Slji[th[t] ] , tfc ' [t] J-w° (- Siii[th[t] ] it[t] + Cos [th[t] ] z [t] ) t- 

2 (l + V (C«M![tJi[t]] if[t] + S±iiLtlk[t]] =[t]) = )*'''] / [i + n*3.[t]^ +^ i[t]^) , 
ie[0] = 0, z[0] = a, tTi[0] = oj , {x, th) , {t, 0, 20}j r 

n - iJcLt] , zLt] } ; 

R = {{Cos [tli[t] ] , SiiiLth[tJ ] } , f- Sin[tli[t] 1 , Cos [th [t ] ] } } ? 
f = Evaluate [ - E . a / . s 2] ; 

pC = Fax awtrlcFlot [f , {t, D, 2«} , PlotKange -> KLl] t 
HajixpiLLatc [Shov[ { Pof ametrxcPlot [ 

Evaluate [Elatten [ft / . a2 / . t -t i , 1] . a • flatten [ - B.- a / . s2 / . t i , 1] / . s2] , 
{t, 0, i} , HotKajuge -» 3111 , PlotStja« - Urective [Red , ThxclcHss [ . 002 ] ] ], 
p£ , Graphics [Point [Elatt^K [fy.t^i]]]}], [i, 0, t}] 



Figure 3.2. This code produces the video suggested in Figure (1.6) 



PazdnetzicPlDtaD j^Eviauate [{x[t] Cos[Sqrt[ (1 - v) / yJ b] + z [tj Sin[ Sqrt [ fl - v) ^ v] a] 
s, -x[t] Sin[Sqrt f [1 - f) / w] h] f z[t] Cos[Sgrt [(1 - v) / v] a] } / . 
HBSolire [jjc ' [t] = Cos[th[t] ] , z'ft] == SiB[th[tJ ] , 

til ' [t] == |-Sqrt[ (1 - -r) / j-Siii[th[t] ] K[t] .CoB[th[t] ] z[l] ) + 
2 (1+ Sqil[(l - w) /r]^ (CoB[tk[t] ] K[t] + SinLth[t]] z[t]l^)*''j / 
(1+Sqrt[(l - Ti /T]'«[t]SSgrt[(l - v) J x]' z[tj'), x[0] = H, 

z[0] == ^ [-ltVlf4M ) , th[l] = •], i», z, th}, {t, 0, t}]] , 

[t, i, !} , {s, «, tt/ Sgi;t[{l -¥)/¥]}, Mesk-f None, 
ExclvsiDiiBStyle -+ {HonCj ted} , 

PlotStyle -» Di£€ctive[Red, Dpacitj [0. t} , SpecaLaritf [Uiite , 20}} j, {H, 

-.24, 

5}. fv, 
.1, .}}] 



Figure 3.3. This code produces the video suggested in Figure (1.7) 
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m|t|;= CleiT [H, w] 

pi = 71 Cos [u] : 

p2 = ( - 1 + Sqrt [1 . 4 M + E Cos [u] "2] Siii[u] ) / (Sqzt [4 + 4w"2a"3 Cos [uj " 2] ) ; 

2 i/l + fl^ ^1^ ^ 1 + {1^ - f2 j 

tl = Sij^liEr f2 + 1 / . ds - 1 ; 

L l.v^ (51=t£!^) J 

2 1 + »^ fl^ + f 2 Jl^ + il^ - f 2 1 

£J = Siii¥lif]>-[-de fl / . da -I 1 ; 

l+»= (£l=f ?2=) ■' 

A = D [ol, u] * Z + D tijj , u] ' 2 / . 

(B-f (2 + 2M'2ir'-4-w'-2t2(Mw'2-l) Sqrt [1+ (1+2 M) w-I+M'-lw'i]) yw*2, 
n ^ Svt [- 1 + K-w" 2 + Sqrt [1 ♦ (1 t 2 MJ w*2 + H--2w*41 ] / (Spt [i] w) ) ^ 
H = f l"-2 + £2 "2 y . {fl H. pi, 42 H. p2 1 / . 

{Bh (2 + 2H*2Tr'-4»w'-2 + 2(Mw*2-l) S^t [1* (1+2M) w'2+M'-2w-^4]) /w'2, 
n - Sqrt [- 1 + Kw'-2 + Sqrt [1+ (1-2M) w-2+M'-2w*41] J (Sqit [2] wl ) : 

2 ^ 1 + 41^ + ^2 41^ -^l fl^ - ^2 j 

;f = da Sart [A y y . (de -. 1 y . {^i ^ pi , ?2 p2 ) y . 

^ l»w^ (Jl^ + f2") J 

{B-t (2 + 2M''2ll''4tw"I + 2(Mw-2-l| S^t [.1+ (1+2M) vi-'-2+M'-2w--4]) /-ur'^, 
n S^t [- 1 + Mir" 2 + iqrt [It (1.2 Ml w"2+M'-2w"4]] y ( S'pt [2] w| } ; 
yO = 0^ 

H ^ tT"" 2; 
ClCAT [w] 

FlniiKoot. [ Nlnt-eqirate [* /. {w -> ir) , {u, , 2 jr) ] - 5 n / 2 = • , {t , »} ] ; 
w = V / . 

Tl = Nlntegiate [Sqrt [A / ;i] , {m, 0. 2 jrf] ; 
so = NBSolve [|>c ' [t] == Cos [tli[t] ] , 

z ' rt] == Sln[tK[t] ] , tk ' [t] == (- Sin[t}i[t] J it[t] + Cos Lth[t] J z [t] 1 + 

2 [1 nr' (C»s [th[t ] J ic[t] 4 Sin[U[t1 ] z[tj ) ^) / (l + w^3i[tl^iw^z[t]^), 

it[0] = e, z[OJ = rO , th[«J == *| , 13C, z, tk) , it, 0, 4 Tl}] .: 

VazaKtricFlDt :Evaliiiit.e [(K[t] , z[.t] } /. so], (t, •, 11} , PlatStyle -» Birective [Ked, Th±ck] ] 
PazaKtricFlot [ {{ Evaluate [{ st [t] , z[t] J y . so] } , {.ETalii^te [ tic[t + Tl] , z [t + Tl] } / . bo] ] , 
(Evaluate [(x[t + 2Tl], z[t.2 Tl] } / . so] ) , (Evalute [ { x [t t 3 Tl] , z[t + 3T1] } y . so] }} , 
{.t , 0, Tl) , (PlotXtf^le -> {Directive [Ked, Tkick] , Directive [Creen, Ihick] , 
Directive [Blue, TKick] , Directive [Yellew, Thick]}}] 
grl = ParBwiricPlotSD [{ (x[t] Cos [wv] <- z [t] Siii[wv], v, -x[t] Sin [w v] -i- z [t] Cos [w v] ) y . so} , 
(t, *, Tl} , [v, *, .7}, Hesh-tNone, Exclusions Strle - {Hone, Ked], 
PlotStyle Directive [Bledj Opacity [ . S ] , Specvlazity [Uhlte , 2D]]j7 

qrJ = Pam^ti^icPlotaD [[{x[t +T1] Cos[wv] + i [t + Tl] Siii.[ wv] , 
T, -x[t +T1] Sim [wv] - z[t + Tl] Cob [wt] } / . so} , 
[t, •, Tl] , [v, •, .1], Hesh -t None, Exclusions St:Fle [Kone, Creen], 
PlotStyle Directive [Green, Opacitj- [0 . t] , Specslaritj^ [White , 20]]]^ 
gr3 = ParawtricPlotaD [{{x[t 1- 2T1] Cob[wt] +z[t + 2Il] Sin[wv], 
V, -x[t + 2 Tl] Sin[wv] + z [t t 2 Tl] Cos [wv] ) y . so] , 
[t, Tl} , {v, a, .1}, Mesh None, ExcluaionsStf le {Hone, Blue}, 
PlotStyle -» Birective [Blue , Opacity [D . t] , Specularity [White , 2P]]]; 
gr4 ^ ParawtricPlotaD [( {x[t 4 3 Tl] Cos [-H v] +z[t + 3Tl] Sin[wv], 
V, -Jt[t + J Tl] Sin[wv] + z [t + 3 Tl] Cos [wv] ] y . so ] , 
{t, •, Tl} , (V, a, .1], Hesh -I Hone, ExclusionsStyle {Hone, Tellair} , 
PlotStyle Directive [Yellow, Opacity [4. t] , Specularity [White , 20}]] r 
Show[{grl, qr2, gr3 , gi-4} , PlotRanqe -» All] 
Clear [M] 



Figure 3.4. This code produces the picture in Figure (1.8) 
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ini[l];= 

Hanipulate ^ 

jpajraKtricPlotan [Evaluate [ {x [t] , Sin [v] z [t.] , Cos [m] z [t] } /. NOSdIve [ 

{Jt ' [t] == Cos [th[ 1 1 1 , z ' [t] == Si»[tli[t] 1 , th ' [t] == 2 + Cos [tli[t] ] y z [t] , x[0] == , e [0] = 
(Sgrt[4 (-Sqrt[v] / (1+ Sg^t[v] ) "2) .1] * 1) / 2 , th[*] == tt), {x, tJi}„ {t, 0, 5)] ] , 
{t, •, 5}, (11, 0, 2 k}, Hesh -• None, ExclusionsSt^i^le -> {Hvne, Ke4}, 
PlotStyle -» Directive [Green , Opacity [0 . tj , Specularity [White , 20]]], 
Pai'ametncVlatSD [Evaluate [[x[t] , Slm[u] z[t] , Cos [uj z [t] } / . HBSalve [ 

{it ' [t] == Cos [th[ 1 1 ] , K ' [t] == Si»[tfc[t] ] , tk ' [fe] ==2 +Cos [th[t] ] / z [t] , x[ll] == , z [0] = 
(S^t; [4 (Sqrt [»] / ( 1 - Sqrt [»-])"-2)+l]+l)/2, tli[0] = jr), (at, z,th), {t, 0, 15}]], 
{t, •, 15), (u, 2 k}, Hesh -I None , Exclusions St^le -> (None, Bed}, 
PlotStyle -t Direetive [ Green , Opacity [ . t ] , Specularity [White , 20]]] 

, Paf awtricPlotai ^Evaluate [ ( x [t ] Cos [Sqrt [(1 - v) /v] s] Sin [ Sgrt t ( 1 - v] / v] s] , 

s, -»[t] Su[Sqrt [(1 - v) / -r] s] + z[t] Cos [Sqrt[(l - v) / v] s]) /. 
NVSolve |^|x ' [t] = Cos [th[t] ] , z ' [t] == Sin[tli[t] ] , 

th ' [t] == J- Sg^t [(1 - V) / ir]^ ( - Sin[th[t] ] x[t] + Cos [th[t.] ] z [t] ) + 
2 [l + Sqrt [ (1 - V) / v] ^ (Cos [th[t] ] ii[t] + Sin[th[t] ] z[t] ) ^) ^"j 
(1 + Sqrt [ ( 1 - V) / V] ^ x[t] ^ + Sqrt [ (1 - t) / t]^ z[t]') , 3c[0] == 0, 
z[9] == 0, tli[0] == o|, {X, z, thf, {t, 0, 8)]] , {t, •, 8}, 
{s, •, ;T/Sgrt[(l - v) / v] ) , Heah -f None , Exclusions Style -r {Hone, Ked) , 
PlotStyle -r Directi-re [Bed, Opacity [D . 8] , Specalaiity [Hkite , 20]]jj, (v, .01, .94}j 

Figure 3.5. This code produces the video suggested in Figure (1.15) 
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